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ANALYSIS OF A VISCOSITY MODEL FOR CONCENTRATED POLYMERS 


MIROSLAV BULICEK, PIOTR GWIAZDA, ENDRE SULI, AND AGNIESZKA SWIERCZEWSKA-GWIAZDA 


Abstract. The paper is concerned with a class of mathematical models for polymeric fluids, which involves the cou¬ 
pling of the Navier-Stokes equations for a viscous, incompressible, constant-density fluid with a parabolic-hyperbolic 
integro-differential equation describing the evolution of the polymer distribution function in the solvent, and a par¬ 
abolic integro-differential equation for the evolution of the monomer density function in the solvent. The viscosity 
coefficient, appearing in the balance of linear momentum equation in the Navier-Stokes system, includes dependence 
on the shear-rate as well as on the weight-averaged polymer chain length. The system of partial differential equations 
under consideration captures the impact of polymerization and depolymerization effects on the viscosity of the fluid. 
We prove the existence of global-in-time, large-data weak solutions under fairly general hypotheses. 


1. Viscosity models for polymers - formulation of the problem 

Contemporary approaches to the modelling of polymeric fluids have exploited multi-scale descriptions in an 
essential way. Mathematical models have thus been built by coupling systems of partial differential equations 
describing the motion of the solvent with equations that track the evolution of the microscopic behaviour of the solute 
in the solution. In this paper we focus on the modelling and the analysis of two phenomena: we wish to explore how 
the rheological properties of the fluid are affected by the presence of, possibly very long, chains of macromolecules; 
and, second, we aim to investigate the possible mechanisms for polymerization and depolymerization effects that 
can also heavily depend on the properties of the flow. The solvent is considered to be an incompressible fluid and 
we assume that the possible changes of the density are negligible relative to other phenomena. Thus, the underlying 
system of equations under consideration consists of the balance of linear momentum and the incompressibility 
constraint, i.e., 

d t (gv) + div x (gv®v) -div^T = gf, 

diva, v = 0. 

These equations are assumed to be satisfied in a time-space cylinder Q := (0 ,T) x 12, with 12 C K d , d £ {2,3}, 
being the physical flow domain. Here, v : Q — > denotes the velocity of the solvent, g is the constant density, 
/ : Q — > represents the density of volume forces and T : Q — > M. dxd is the Cauchy stress. In order to close the 
system CD, one also needs to relate the Cauchy stress to other quantities describing the qualitative behaviour of 
the material. Despite their importance, in the present paper we shall, for the sake of simplicity, neglect all thermal 
effects and will focus instead on mechanical properties of the fluid in the isothermal setting. In particular, we 
are interested in models that link the Cauchy stress to the shear rate and to the contribution of the microscopic 
polymer molecules, in order to investigate how the level of interaction between polymer chains, their lengths, and 
their configuration influence the Cauchy stress. 

That the presence of macromolecules in a solvent dramatically changes the properties of the flow was already 
observed in [13] . and has been, ever since, the focus of mathematical models. In general terms, one can find two 
different approaches to the problem: the empirical (rheological) one, where one typically fits the observed data in 
the model and introduces empirical formulae for the viscosity, u, based on the ratio of the magnitude of the Cauchy 
stress and the magnitude of the shear rate; and a second one, based on modelling the microstructure of the polymer 
under consideration using tools from statistical mechanics. The second approach has been particularly common in 
the case of dilute polymers, i.e., when the concentration of polymer molecules in the solvent is very low, so that 
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the interactions of the macromolecules can be neglected. In such models the presence of macromolecules in the 
solvent is typically modelled by additively supplementing the Cauchy stress tensor with an elastic stress tensor. 
This then leads to a closed system of partial differential equations; for a derivation of the model and further related 
mathematical results and references we refer the reader to Barrett & Siili Euan mm and Bulicek, Malek & Siili 
[9]; see also pT8] and [I]. 

On the other hand, in the case of concentrated polymers one has to take into account the interactions of the 
polymer molecules (see, for example, |12]). Here we shall be concerned with models in which the only important 
property of the microstructure that influences the macroscopic flow equations is the polymer chain length, in the 
sense that the viscosity coefficient appearing in the balance of linear momentum equation is considered to be a 
function of the polymer chain length, resulting in a family of viscosity models. There is an evident contrast between 
mathematical studies of models that ignore the dependence of the viscosity coefficient on the lengths of polymer 
chains, and the vast set of experimental data indicating the influence of the polymer chain length on the viscosity, 
see m ini m da on among others. In particular, none of the approaches mentioned above considered the impact of 
polymerization and depolymerization effects on the flow properties. The goal of this paper is therefore to explore a 
class of viscosity models that incorporate the influence of polymerization and depolymerization on the macroscopic 
flow variables. Specifically, we provide a rigorous proof of the existence of a global in time, large-data weak solution 
to the class of models under consideration. 


1.1. Viscosity models: the relationship between the viscosity, the shear rate and the polymer chain 
length. It was experimentally observed more than seventy years ago already (cf. m, for example,) that the 
properties of the solvent heavily, and nonlinearly, depend on both the number of polymer chains in the solvent 
and the polymer chain lengths. Flory had studied the relationship between the viscosity of the solvent and the 
number of polymer chains in the solvent for linear polyesters at constant temperature, pressure and shear-rate. 
Guided by experimental evidence he proposed that for linear polyesters the logarithm of the viscosity should be 
a linear function of the square-root of the so-called weight-averaged chain length. Moreover, it was also observed 
that this relationship is independent of the type of distribution function for the species in the polymer. To be more 
concrete, if we denote by ip(r) the distribution function of polymer chains of length r (meaning that the polymer 
chain consists of r monomers) and denoting the minimal length of a polymer by ro, the total weight of the polymer 
can be considered to be proportional to 


rip{r) dr. 


This is indeed the case if one considers ro 1, but must be corrected for lower values of ro, see |13j . Thus, the 
weight-averaged chain length can be defined as 




IrT r2 ^( r ) dr 

dr ' 


It was empirically shown in m that 



in the case when one considers linear polyesters; see [14J for a more detailed account. These experiments indicate 
that, generally, one does not need to consider specific values of the distribution function but that its average in the 
above sense will suffice from the point of view of accurate macroscopic modelling. This was also conjectured to be 
the case in many other situations, and quite recently also in the case of amines at high pressures, see [20] . In that 
paper, the authors showed that the logarithm of the viscosity is not necessarily related to the square-root of the 
averaged chain length but it can be a more general function, though still dependent only on the averaged properties 
and not on the fine details of 

Further, it was common belief (based on laboratory experiments at nonconstant shear rates) that many such 
phenomena can be explained by using a shear-dependent viscosity instead of a polymer-length-dependent viscosity. 
However, this is not the case for concentrated polymers. For example, in mi, the authors studied the simultaneous 
influence of the shear-rate and the polymer-length on the viscosity for polydisperse polymer melts. It was observed 
that, up to a certain critical value of the shear rate, the fluid behaves as a Newtonian fluid with viscosity depending 
only on the chain length, while for large shear-rates significant shear-tliinning was observed. Moreover, it was also 
shown that, asymptotically, the viscosity depends only on the shear-rate and the influence of the chain length in 
large shear-rate regions can be neglected. In addition, it also transpired that the weight-averaged chain length was 
not an appropriate quantity in the cases considered, and the authors suggested instead the following general form 
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of weighting: 

( 2 ) 


*P := 


IrT w(r)np(r) dr 
f r °° rip(r) dr 


where the function ui represents the most significant lengths of the polymer chains. Surprisingly, during the last 
decade, it was observed that there exist materials (e.g. polystyrene-decalin, polyethylcn) that exhibit the opposite 
phenomenon: shear-thickening, see m ■ For all these reasons, we consider a much more general form of the Cauchy 
stress: 

T := S — 5 1, 

where q is the pressure and S represents the viscous part of the Cauchy stress; we further assume that S is of the 
form 


(3) 


S(-0, D^n) := v(ip, |Da;u|) D x v , 


where is the symmetric velocity gradient and ip is & general weight-averaged chain length function of the form 
©, which one can simplify (without loss of generality) to the form 

nOO 

ip := l(r)ip{r) dr. 

Jr 0 

We shall not specify the particular form of v, but will instead admit a general class of viscosities so as to enable 
the consideration of both shear-thinning and shear-thickening fluids. 

1.2. Polymerization models. In the above system the microscopic parameter that is taken into account is the 
polymer chain length, and we describe the process of elongation of these chains (by binding free monomers to 
polymer chains) as well as the process of breaking longer chains into shorter ones. We shall view this microscopic 
structure as an evolution of two populations: the population of polymer chains and the population of monomers. 
This kind of description is akin to the model for prion dynamics considered in HU; see also mm- As the model 
formulated in the current paper can be seen as an extension of the prion proliferation model in those papers to 
the case with spatial effects (transport by a solvent and spatial diffusion) we shall discuss it in more detail. The 
authors of mmm study how the healthy prion protein and infectious prion protein populations interact in an 
infected organism. The infectious prions are abnormal pathogenic conformations of the normal ones. The proposed 
model considers the infectious prion proteins to be a polymeric form of normal prion proteins. Polymers of infected 
prions can split into shorter chains. Such a splitting transforms one infectious polymer into two shorter infectious 
polymers, which can then attach again to normal prion proteins. However, when the length of a part of a split 
polymer falls below a certain critical value, it immediately degrades into a normal prion protein, namely a monomer. 
By ip we denote the distribution function of polymer chains of length r > ro that satisfy the following equation: 


(4) 


d t ip(t,r) +T(p(t)d r ip(t,r) = -/3(r)ip(t,r) 


/3(f)«:(r, f)ip(t , r) dr. 


The term T(p(t)d r ip(t, r) represents the gain in length of polymer chains due to polymerization with rate r > 0, /3(r) 
is the fragmentation rate, namely the length-dependent likelihood of splitting of polymers to monomers, n{r,r) is 
the probability that a polymer chain will split into two shorter polymer chains of length r and f — r, respectively, 
the term — f3(r)ip(t , r) is the loss of polymer chains subject to the splitting rate /3(r), and the last term is the count 
of all the polymer chains of length r resulting from the splitting of polymer chains of length greater than r. 

The evolution of monomers is described by means of the function <p(t), which is the concentration of free monomers 
at time t, satisfying the equation 


(5) 


pr o poo poo 

—<j>(t) = 2 / r P(r)K(r,f)ip(t,r)drdr — (p(t) / Tip(t,r)dr. 
dt Jo Jr Jr„ 


The first term on the right-hand side represents the monomers gained when a polymer chain splits with at least 
one polymer chain shorter than the minimum length ro, while the second term is the loss of monomers as they are 
polymerized. 

1.3. The complete model. We conclude this introductory section with the precise statement of the complete 
system of equations under consideration. For a given Lipschitz domain 12 C M d , d €E {2, 3}, and a given final time 
T > 0, we consider the balance of linear momentum and the incompressibility constraint on Q := (0 ,T) x 12 in the 
form (after scaling by the constant density g) 

^ dtv(t, x ) + di v x (v(t, x) <g> v(t, x)) + V x q(t, x) - div x S(^(2, x),D x v(t, x)) = /, 

di v x v(t, x) = 0, 
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where v : Q —>• R d denotes the velocity of the fluid (solvent), q : Q —> R is the pressure, and / : Q —> R d is the 
density of the external body forces. The viscous part of the Cauchy stress S : Q —> R dxd is given by the formula 

(7) S (0(2, x), D x v(t, x)) := v(ip{t, x), |D x v(t, x)|) D x v{t, x). 

Here. D„z; denotes the symmetric velocity gradient, i.e., D x v := |(V x t; + (V a; , i;) T ), and the viscosity coefficient 
v : R + x R + —> R + is allowed to depend on the shear-rate (D^ul and on the averaged polymer distribution function 
0 : Q —> R + , defined by 

(8) ip(t,x):= 'y(r)ip(t,x,r)dr, 

J Mq 


where Ro := (7-0,00) and 7 : Ro —> R+ is a continuous nonnegative function, representing the weight function 
associated with the averaging of polymer chain lengths. The distribution function ip : Q x Ro —>■ R+, which in the 
absence of fluid motion satisfies (ED i is now assumed to satisfy the following equation: 


(9) 


<9 t 0(t, x, r) + v(t, x) • V x ip(t, x, r) + r(r)0(f, x)d r ip(t, x, r) - A(r)A x ip(t, x, r) 


pOO 

= - /3(r, v , D x v)ip(t, x, r) + 2 / /3(f, v, D x v)n(r, r)%p(t, x, f) df, 


in Q x Ro- There are therefore two additional terms compared with the convective/transport term v ■ X7 x i/j due 
to the motion of the solvent, and the diffusion term AA x i/j associated with the Brownian force acting on the polymer 
molecules immersed in the solvent. The parameter ro G [0, 00) is a fixed minimal length of a polymer molecule, 
A{r) : Ro —> R+ denotes the rate of diffusion for a particular value of r and is assumed here to be a nonincreasing 
function of r, t : [ro,oo) —> R+ is the polymerization rate, /3 : Ro x R d x S, dxd — > R + is the fragmentation rate 
of polymer chains, which, unlike (U), will also be allowed to depend on the macroscopic quantities in the model, 
namely on the fluid velocity and the shear rate, and, finally, n(r, r) denotes the probability that a polymer molecule 
of length f will split into two polymer molecules of lengths r and f — r, respectively. The function <j> appearing in 
© is the concentration of free monomers satisfying, in the absence of fluid motion, the identity ©, and which, in 
the presence of a moving solvent, is assumed to satisfy the following equation: 


d t 4>{t, x) + v(t, x) • V^f, x) - A 0 A x (j>(t, x) 

(10) f°° f ro f°° 

= — (f>(t,x) / d r (rT(r))i/j(t, x, r) dr + 2 / r (3(r,v,D x v)k{t, x, f ) df dr 

J ro J 0 J ro 


in the space time cylinder Q. Here, the rate of diffusion Aq (caused by Brownian noise) is assumed to be a positive 
constant and the additional transport term v ■ V x <j> is due to the flow of the solvent. 

We shall further assume that the velocity satisfies a Navier slip boundary condition, i.e., 


( 11 ) 


v ■ n = 0 on (9S2, 

(S n) T = — a*v on <5>ST, 


where a* > 0 is the domain-wall friction coefficient, n is the unit outward normal vector to dfl, and for any v we 
have denoted by v T := v — (v ■ n)n the projection of v on the tangent hyperplane to the boundary. The function 
c t> is assumed to satisfy a homogeneous Neumann boundary condition with respect to the x variable, i.e., 


( 12 ) 


V x <f> ■ n = 0 on (912, 


and for ip we also prescribe a homogeneous Neumann boundary condition with respect to the x variable, i.e., 

(13) V x ip ■ n = 0 on dtt. 


We shall further assume that ip vanishes at infinity with respect to r, i.e., 

(14) lim ip(t,x,r) = 0. 

' ' r—too 


Finally, to complete the statement of the problem, we prescribe the following set of initial conditions: for the 
velocity field we assume that 


(15) t>(0, x) = Vq{x) in 12, div x v 0 = 0 in 12, 
and for <p and ip we assume that 

(16) </>(0, x) = <po(x) in 12, 

(17) 0(0, x, r) = ipo(x,r) in 12 x (r 0 ,oo), 


Vq ■ n = 0 on <912, 


00 > 0 , 
ipo > 0. 
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1.4. Assumptions on the parameters. In what follows K will signify a universal positive constant. We adopt 
the following assumptions on the data: 

(Al) The diffusion rate A : lo -> R + is a continuous nonincreasing strictly positive function defined on Ro such 
that huL^oo A(r) = 0. 

(A2) The polymerization rate r is a nondecreasing bounded globally Lipschitz continuous function such that 
T ( r o) = 0, r'(r 0 ) > 0, and there exists a positive constant K such that 

( 18 ) K~ 1 ro < r(r) + r T r (r) and r(r) + /(r) + rr'(r) + — < £ 

r 

(A3) The fragmentation rate /? : Ro x R d x M. dxd —> R + is a smooth bounded function, which is increasing with 
respect to its first variable, and satisfies, for all (r, v , D) £ R 0 x R d x R dxd , 

( 19 ) 0 < /3(r,v, D) < K. 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 


Moreover, we require that the function 77 , defined by 


d r /?(r,u, D) 

77 (f) := sup — -—, 

ueitDeR d xd p{r,u,u) 


is measurable and nonnegative, and 

f‘00 

0 < (1 + r)r}{r) < K, / 77 (f) df < K. 

J r n 


f 


(A4) The kernel function k is, for simplicity, defined as 


. . — if f > f 0 and 0 < r < r, 

n(r, r) := \ r 

0 otherwise. 


It therefore follows that, for all a > 0, 


r a 1 n(r, r)dr = 


for r > f 0 . 


(A5) There exists a positive real number 9 > 0 such that, for all r £ [ro, 00 ), one has 
(24) 


7 (f) < K(1 + r) s . 

that v : R + x R + —>• R + : 
i~i£ R dxd fulfilling £ 7 A £, one has (cf. ©) 


(A 6 ) We assume that v : R + x R + —> R+ is a continuous function such that, for some p > and for all 


(25) 


\S(;Z)\<K(l + \Z\y-\ 
SM)-£> K-'W-K, 
(S(*, £) — S(-, £)) •(£ — £)> 0. 


1.5. The main result. Before introducing the definition of weak solution to problem © (ED we shall summarize 
our notational conventions. We denote by T £ (0,oo) the length of the time interval and by fi C R d , d £ {2,3}, 
a bounded domain in R d with C 1,:L -boundary dil: we then write £ C 1 ’ 1 . We also set Q := (0,T) x f l and 
T := (0, T) x dfl. 

For p £ [l,oo] we define the Lebesgue space L p ( fl) and the Sobolev spaces W 1,p (Jl) in the usual way, and we 
denote the trace on dfl of a Sobolev function u by tr u. If A is a Banach space, then X d := X x • • • x X and 
we write X* for the dual space of A; L p (0,T; A) denotes the Bochner space of A-valued L p functions defined on 
(0,T). For (scalar, vector- or tensor-valued) functions / and g and • signifying the product of real numbers, scalar 
product of vectors, or scalar product of tensors, as the case may be, we shall write 

{f,9)-= [ f(x)-g(x)dx if f-g£L 1 (n), 

J n 

(f,9)on~ f f(S)-g(S)dS ii f ■ g £ L\dft), 

JdQ. 

{gJ) '■= (gj)x*,x if / G A and g £ X*. 

We also require the space C wea k(0,T; L p (tt)) consisting of all u £ L°°(0, T; L p (fl)) such that (u(t),<p) £ C([0,T]) for 
all ip £ C(H). 


$LaTeX: 2016/1/8$ 





6 


M. BULICEK, P. GWIAZDA, E. SULI, AND A. SWIERCZEWSKA-GWIAZDA 


We introduce certain subspaces (and their dual spaces) of the space of vector-valued Sobolev functions from 
fV 1,p (12) d , which have zero normal component on the boundary. First, we define, in the usual way, for any p £ [ 1, oo), 

LP n,div := i v e di Va: t> = 0}" " P , 

where by 22(12) we mean the set of all infinitely smooth functions with a compact support in the set 12. Then, by 
V and Vdiv we denote 


V := {v £ W d+2 ' 2 (VL) d ] v ■ n = 0 on <912}, Vdiv := V fl L 2 ndiv . 

Note that V C and therefore we can finally introduce the following spaces for any p £ [l,oo), and 

p' :=p/{p- 1): 

W^:=V IMkp , W .:= {W^Y, 

Wjdiv W~X := (W&)'. 

Moreover, for any a > 1 we introduce the space LP, (R 0 ) : = {v G T p (Ro) '■ f Ro r a p p (r ) dr < oo}; analogously, we let 
L p (fl x M 0 ) := {p £ L p (12 x R 0 ) : f Q f Ro r a p p (x, r) dr dx < oo}, R 0 := (r 0 ,oo). 

Theorem 1.1. Let the assumptions (A1)-(A6) be satisfied. Then, for any 12 G C 1 ' 1 , T G (0,oo), any v d , f, and 
any nonnegative to, to satisfying 

(26) v 0 GL^ V , /Gf7'(0,T;^'), to G ^(12; L\. (M 0 )) U L 2 (12; L§(Ro)), to G L°°(n), 
with some 0} > 6 and Of > 1, there exists a quadruple (q,v,t,t) an d Q* > 1 such that 

(27) «G2A*(Q), 

(28) v G C wea k{ 0, T 1 ; L n ^ j v )ni p ( 0 ,T;<; p div ), 

(29) ^ G L°°(0, T; L§(f2 x R 0 )) n L 2 (0, T; W 1 ' 2 ^))) n ^(O, T; T 1 (f2; L\, (M 0 ))), 

(30) G i“((0,T) x 12) ni 2 (0,T; W 1,2 (fi)), 

with 8 t v G L«’(0,T;W"^f), d t ip G L 9 *(0,T; (W 1,2 (]R 0 x 12) n x 12))*) and d t t G T 2 (0,T; W 1 ’ 2 (12)), 

which attains the initial conditions in the following sense: 

(31) lim \\v(t) - v 0 \\l + || tit) - to IIi + Wtif) ~ to II 2 = 0, 

t— HJ-(- 

and solves, for almost all t G (0,T), the equations ([HI) (1131) in the following sense: for all w G W^ 1 such that 
\/ x w G L°°{fl) dxd , 

(32) (d t v, w) + (S, V x w) - (v <g> v, V x w) + a*(v, w) an = (/, w) + (q, div x w ); 


for all p G W 1 ’°°(12;22(R)), 


(33) 


(dtt, t) - ( Vt , VxP) - (d r (Tip), tt) + (M r )^xt, V x t ) 

= ~{Pt,v) + 2 i^J ^Ktdr,pj; 


and for all z G W 1,2 (12), 


(34) 


z) - (vt, V^z) + A 0 (V x t, V x z) 



d r (rr)t dr, z 


\ 

/ r r 0 

r-OO 

)+2 

r 

/ 

\J 0 

Jr 0 


(3nt dr dr, 



A key difficulty in the mathematical analysis of the problem under consideration is that the partial differential 
equation (|9|). governing the evolution of the distribution of polymer chains, is a hyperbolic equation with respect 
to the variable r with a nonlocal term, which is nonlinearly coupled to the evolution equations © and (flT)l) . 
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2. Analytical framework 


An essential part of the existence proof in the fluid part of the problem relies on Lipschitz approximations of 
Bochner functions taking values in Sobolev spaces. We recall from [8] the following lemma, which collects the 
properties of the approximation in a simplified setting (omitting the generality of Orlicz spaces used in [8]). 


Lemma 2.1. Let f2 C d £ {2, 3}, be a bounded open set, let T > 0 be the length of the time interval and suppose 

notion H and arbitrary seqv 

a n := |H"| + IHI and b n := \D x u r 


that i + y = 1, with p £ (1, oo). For any function H and arbitrary sequences {u n }ff = i and {H n }^L 1; we consider 


such that, for a certain C* > 1, 

[ \a n \ p ' + \b n \ p dxdt + ess.sup te(0iT) ||u n (f)||^ < C*, 

(35) Jq 

u n —> 0 a.e. in Q := (0, T) x Q. 

In addition, let {G"}~ ! and {/”} r T =1 be such that G" is symmetric and 

(36) G n 0 strongly in L 1 (Q ) dxd , 

(37) f n -*0 strongly in L 1 (Q) d , 

and suppose that the following identity holds in T>'(Q) d : 

(38) 8 t u n + divz(H” - H + G n ) = 

Then, there exists a /3 > 0 such that, for arbitrary Q k CC Q and for arbitrary X* £ (pp ^ 1 , oo) and arbitrary k £ N, 
there exists a sequence of {X^}^ =1 , a sequence of open sets C Q, and a sequence bounded 

in L < ^ c (0,T-,W^^°(Il) d ), such that, for any 1 < s < oo, 


(39) 

(40) 

(41) 

(42) 

(43) 

Moreover, 

(44) 


u 


K € 

n,k _. 


p -i k 

X*,p-—(X*) P 


\D x (u n ’ k )\\ L ~ {Qh) <c(h,n)x n k , 


U ’ = u 


Vn £ N, 

strongly in L s (Qh ) d , 
in Q h \ E%, 


limsup | Qh (~l E k \ < C(h, f2) 


C* 


(A*)? 


limsup f (|H"| + |H|) \D x (u n ' k )\dxdt < C{h,C*) 

n—too jQ h OE™ 


1 


(A*) 


,p -1 


1 


and the following bound holds for all g £ T>{Qh): 

( 45 ) - liminf j\dtu n , u n ’ k g) d t < C(g, h, C*) + \ 

3. Uniform a priori estimates 

In this section, we shall formally derive the uniform estimates that will play a crucial role in the proof of existence 
of weak solutions to the problem. All of the statements below are valid for sufficiently smooth solutions. First, 
we recall a minimum principle. However, we will prove it for the following slightly modified problem (we shall 
only indicate the dependence of functions on the variable r, except in cases when it is necessary to emphasize the 
dependence on the other independent variables as well): 

/»00 

(46) d t if + v ■ V x if — AA x if = — Pip — T(j>d r if + + 2 / f3(f, -)K,(r, r)ip+(-, f) dr, 

J r 

r-OO /»ro /*oo 

( 47 ) d t <\> + v ■ V x (f> - A q A x 4> = -(f / d r (rT)ip + dr + 2 / r /3(r,-)K(r,r)i/> + (-,r)drdr, 

J 7*0 «/ 0 J Tq 

where := max(i/>,0). 

Lemma 3.1 (Minimum principle for if and </>). Let ifo and </>o be nonnegative; then, the functions if and <f that 
solve the coupled system m and 63 ) are also nonnegative. 
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It directly follows from Lemma 13.11 that we can replace ip + by ip and then (HSl) , (1771) reduce to m, m- 


Proof. We test the equation (1461) by ip-, where if- = min('0,O) and integrate with respect to x and r. Using the 
Neumann boundary condition on -0 we get, after integration by parts, that 


(48) 


+ J n J r ( A (')I V V’-| 2 +/3(uOIV>-| 2 ) drdx 


/» /*00 /» n OO n OO 

- / / (diva, v)\ip-\ 2 + T(pip-d r ip + dr dx + 2 / / ip- f3(f, ■)n(r,f)ip + (-,f) dr dr dx 
^ JQ Jro JO, J rn Jr 


Q Jr 0 


where we have used the fact that v ■ n = 0 on dfl. Then, the first term on the right-hand side is identically zero 
since div^, v = 0. The same is true for the second term on the right-hand side since ip-d r ip+ = 0. Finally, since /? 
is nonnegative, we see that the last term is nonpositive. Consequently, we get 

and since we have assumed that ipo > 0, we deduce that ip > 0 almost everywhere. Using a similar procedure, we 
obtain an identical result also for cp. □ 

Our next estimate is a maximum principle for cp. 

Lemma 3.2 (Maximum principle for cp ). There exists a constant C, depending only on K, such that , if <po € 
then 


(49) 


ess.sup te(0T) 


< max (A' 2 ,110o 11 oo)- 


Proof. We begin the proof with a pointwise bound on the right-hand side of m- Using the nonnegativity of ip 
and (p (Lemma 13.II) and the assumptions on r and (3 and the definition of k, we observe that 


(•OO 


— (p d r (rT)ipdr + 2 


/3(f, -)«;(r, r)ip(r) df dr 


Jr 0 •> U •> To 

poo pr 

< —<p / (t + rd r r)ip d r + K 

Jr 0 JO Jr 0 

r°° r r o r°° j,(f 

< —K~ 1 r 0 (p / ipdr + K r - : 

J ro JO .//(j P> 

Kr 0 


ip(f) 


df dr 


r 

■ df dr 


p OO 

= —K~ 1 rQCp / ip dr 
J r 0 


ip dr 


' r O 


< -K~ x r o 


< -AT V 0 


/ ip{r) dr J {(p — A' 2 ) 
J r Q / 

J ip(r) dr^ ((p — M ), 


where M is defined as M := max(A' 2 ,110o 11 oo) - Hence, by multiplying (H71) with (cp — M) + and integrating over U, 
we get (using integration by parts, the Neumann data and the fact that the velocity is a solenoidal function) that 


-||(^-M) + || 2 <0. 


Consequently, we immediately arrive at OH). 

The next result concerns conservation of mass. 


□ 


Lemma 3.3 (Conservation of mass). Let the pair of functions (ip, (p) be a solution to dHb (1101) : then, the following 
identity holds: 

d 


(50) 

Consequently, for a.a. t £ (0 ,T), 


dt 


(p(t,x)dx+ / r 
! J ro 


ip(t, x, r) dx dr 


= 0. 


p p OO p p p OO p 

(51) / (p(t,x)dx+ / r ip(t, x, r) dxdr = / cpo(x)dx+ r ipo(x, r) dx dr =: Eq. 

JQ J ro J Q J Q J ro J Q 
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Proof. First, we integrate © over Cl. Since div x v = 0, the transport term is equal to div x (vip) and then both the 
transport and diffusion terms vanish thanks to v ■ n = 0 on dCl and the zero Neumann boundary condition for ip. 
Hence, we obtain 


(52) 


dt / ip(t, x, r) dx + / <p(t, x)Td r ip{t, x, r) da; 

J Q J £2 

= — / /?(r, -)ip{t, x, r) dx + 2 / / n,{r,f)j3{f^)ip{t,x,f) dr dx. 

J Q JQ Jr 


(53) 


In the next step we multiply (15211 by r and integrate over the interval (ro, oo). Consequently, 
d f°° f f°° f 

— r ip(t, x, r) dx dr + / / rr(p(t,x)d r ip(t,x,r)dxdr 

dt J ro J a Jr 0 J fi 

p OO p pOO p pOO 

= — / r(3(r, -)ip(t, x, r) dx dr + 2 / / r / K(r,f)(3(r,-)ip(t,x,f)drdxdr. 

J ro J Q J ro J Q, J r 

Finally, assuming that ip vanishes sufficiently quickly at infinitj0, we can integrate by parts with respect to r (note 
that since t(tq) = 0 the second boundary term also vanishes); this yields the identity 

d f°° f r°° [ 

— / r I ip(t, x, r) dx dr — / / d r {rr)(p{t,x)ip{t,x,r) dx dr 

dt Jr 0 Jn Jro J a 

p OO p p OO 

/3(r,-)ip(t,x,r)dxdr + 2 / / r / K(r,r)/3(r,-)ip(t,x,f)drdxdr. 

J ro «/ Q J r 


(54) 


' r o Jf2 


Next we integrate m over Cl. Similarly as above, after integration by parts, the transport and the diffusion terms 
vanish and we get 

d 


(55) 


dt 


/ </>(£, x) dx = — / / <j>(t, x)d r (rr)'i/j(t, x, r) dr dx 

J Q J Q J ro 

r r^o noo 

+ 2 / / r f3(f,-)n(r,f)ip(t,x,r)dfdrdx. 

J Q J 0 J ro 


Moreover, using the fact that n(r, f) = 0 for r < ro, we can rewrite the last term to obtain the identity 

/ </>(f, x)dx = — 

t Q J ro 


(56) 


/ <p{t,x)dx = — / / <p(t, x)d r (rT)ip(t, x, r) dr dx 

J Q J Q, J ro 

r r r o /*°° 

+ 2 / / r (3(f,-)n(r,r)ip(t,x,f)drdrdx. 

J f2 J 0 Jr 


Thus, by summing (|54| and (l56l) and using Fubini’s theorem, we have that 
■ I / <p(t,x) dx + 

la Jr 0 


(57) 


dt 


r r r°° > 

/ cp(t, x) dx + / / rip(t, x, r) dr dx 

J Q J Q, Jro J 

p poo p poo poo 

= — / / r/3(r, -)ip{t, x, r) dr da: + 2 / / r K(r,r)/3(f,-)ip(t,x,r)drdrdx. 

J Q J ro J Q Jo Jr 


'Q Jr 0 

Finally, we evaluate the last term. Changing the order of integration and using the fact that k(t, f) = 0 for r < r 0 , 
we deduce that 


/» /* OO /»00 

/ / r re(r, r)/3(r, -)ip(t, x, f) dfdr dx 

J Q, J 0 Jr 

p poo poo 

/ / / Xr>rrn(r,r)P(r, -)ip{t, x, f)dfdrdx 

Jn Jo Jo 

/* /*oo pr 

/ / / rn{r, r)/3(f, -)ip(t, x, f) dr dr dx 

Jn Jo Jo 

n l*oo i nr \ ([ 23 ]) f 

/ / / T/t(r, f) dr I /3(r, x, f) dfda; = / / f/3(f, ^-^(t, x, r) dr dx. 

J Q J ro \J 0 J J Q, J ro 


in JO 
= 2 

= 2 


= 2 

Consequently, we see that the right-hand side of (l57l) is identically zero, and we deduce (l50l) . 
We shall now develop further bounds on the function <j>. 


□ 


1 At the level of these formal computations, it is assumed that i/j vanishes sufficiently rapidly at infinity; the argument will be made 
rigorous in Section [4] by fixing the function space in which the tripe ,</>), understood as a weak solution to the problem, is sought. 
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Lemma 3.4 (Parabolic regularity of ft). Let ipo, ft o be nonnegative, fto 6 L°°(Ll) and E 0 < oo (cf. (15T1) ): then, 


( 60 ) 


(58) / \Vft\ 2 dxdt<C(\\fto\\oo,E 0 ,K,T). 

JQ 

Proof. Let us rewrite equation m as follows: 

(59) dtft + v- V x ft - A 0 Aft = F, 

where F is the right-hand side of ( 11011 . First, we deduce a pointwise bound on F. Using the minimum principle for 
■ft in conjunction with the properties of /3 and r, we get 

/•oo nr o n oo 

F(t, x) = — ft(t, x) / d r (rr)ft(t, x, r) dr + 2 / r j3{r,-)n{r,r)ft(t,x,r) dr dr 
Jr o J 0 J ro 

pro poo pro poo 

<2 / r fi(f,-)n(r,r)'ft(t,x,r) dr dr <2K / r n{r, r)ijj(t, x, r) dr dr 
J 0 J ro J 0 J ro 

poo / pro \ p oo 

= 2K / I / r«(r, r) dr I ftft, x, f) dr < K / roij}(t,x,r) df 

J ro \J 0 J J ro 

poo 

<K r^(t,x,r)dr. 

J ro 

Thus, upon multiplying (15U1) by ft, integrating over f2, using the nonnegativity of ft, partial integration (noting that 
all boundary terms again vanish) together with ( 1511 ) and Lemma 13.21 give: 

1 f ] r r°° 

2 -fa\\ft ( t )\\2 +A 0 \\Vft(t)\\l < F J ft{t,x) J rip(t, x,r) dr dx 

< K\\ft{t)\\oo [ j rft(t,x,r)drdx 

J 17 J ro 

mm 

< KE 0 max(iF 2 , H^olloo)- 

Hence (l58l) directly follows. We note in passing that we have not explicitly indicated the dependence of the constant 
C'GIV'ollcx^ Eq, K, T) appearing in (I581l on A 0 and |U|. □ 

The next step is to establish a bound on high-order moments of the function if. 

Lemma 3.5 (High-order moments of ft). Let ft be nonnegative and suppose that it satisfies <[9|); then, for all a > 0, 


( 61 ) 


) p poo 

ess. sup te(0 ,T) / / r a ft(t,x,r)drdx < C(\\fto\\oo,ct,K,T) / r a ft 0 (x,r) dr dx. 

J 17 J ro J Q J ro 


Proof. We first integrate (0 over H to get 

dt / ft(t,x,r)dx = — / fi(r, -)ft(t, x, r) dx + 2 / / /3(f, ■ )n(r, r)ft(t, x,f) dr dx 

J n J n Jn Jr 

— / r{r)ft{t, x)d r ft(t, x, r) dec. 

Jn 

We then multiply the result by r“, where a > 0, and integrate over (ro, oo) to deduce that 

^ f f r a ft(t,x,r) dr dx = — f f r a f3{r,-)ft(t,x,r) dr dx 
J Qj J ro J Q J ro 

n f OO n OO n n OO 

+ 2 / / r“ / /3(r,-)n(r,r)ft(t,x,r)drdrdx — / / T(r)r a ft{t,x)d r ft(t,x,r) dr dx. 

J Q J ro J r J Q J ro 


dt 
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Next, we evaluate the second term on the right-hand side. Using Fubini’s theorem, we arrive at 

pOO poo poo poo 

2 / r a /3(r, -)n(r, r)if(t, x, f) dr dr = 2 / / r a Xr<rP{r,-)K(r,r)il)(t,x,r)drdr 

J rn J r J ro J ro 


/r 0 Jr 

poo pr 

= 2 / / r a /3(r, r)^(t, x, r) dr dr 

Jr 0 J 7*0 


= 2 

f r 0 
J23l 2 


/ oo / rr \ roo / rr 0 \ 

/ r a ft(r, f) dr I/3(f, x, f) df — 2 / ( / r a «(r, r) dr ) /3(f, x, f) df 


o^a+l /*°° 
0 


a 


- / r a /3(r,-)ij}{t,x,r)dr -— 

1 a + 1 


/3 (r,-)ip(t,x,r) 


dr. 


Hence, we have that 

d ' '■°° 


dt 


r a i/.’(t, x, r) dr dx H- —- 

a + 1 


Q Jro 
i „ r poo 

1 — a ‘ ' 

a + 1 Jq j ro 


P(r, -)i>{t,x,r) 


dr dx 


fi Jro 


t [ [ r a /3(r, -)ip(t, x, r) dr dx — f [ T{r)r a (j)[t, x)d r ip(t, x, r) dr dx. 

1 Jn Jro Jn Jro 


Using the nonnegativity of if, the boundedness of <f> and integration by parts in the last integral, 


1 p poo 1 _ p poo 

— r a i/j(t, x, r) dr dx < - / / r a /3(r, d'tpft, x, r) dr dx 

d t Jn Jro a + 1 7o Jr o 

+ / I d r {r(r)r a )(f>(t, x)ip(t, x, r) dr dx 
J Q J ro 

< C(a, K)(l + ||<^(t)||oo) f f r a, 0(t, x, r) dr dx, 

J Q J ro 

where we have used the assumptions on r and fi. Thus, (1611) follows by using Gronwall’s lemma. □ 

Now we can prove the desired bounds on if). 

Lemma 3.6 (Parabolic-hyperbolic estimates for i/;). Let if solve ©■ Then, for all p £ [3,oo), the following 
inequality holds: 


(62) 


p poo p poo 

.sup te ( 0iT ) / / r p ip 2 {t, x, r) dr dx + / / r p A(r)\S7 x ip{t, x, r)| 2 dr dx dt 

J Q J r o J Q J r q 

< C(p,K, ||</> 0 ||oo) f f r p tfl(x,r)drdx. 

J Q J vq 


Proof. Let a(r) be an arbitrary nonnegative function. We multiply © by a{r)if(t, x, r) and integrate over f1 x 
(ro, oo) to deduce that (note that the term containing v again vanishes) 


Id f ,0 ° 
2 dt 


( 63 ) 


p poo p poo 

/ / a(r)ip 2 (t, x, r) dr dx + / / a(r)H(r)|V x ^(f, x, r )| 2 dr dx 

J Q, J ro J 17 J ro 

/ / a(r)/3(r, -)ip 2 (t, x, r) dr dx 

J Qj J tq 

p poo poo 

2 / / / /3(f, -)/c(r, r)if(t, x, f) dr a(r)if(t, x, r) dr dx 

JQ J ro Jr 

— f f 4>(t, x)T(r)a(r)d r ip 2 {t, x, r) dr dx. 

2 7n Jr 0 
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We begin by focusing on the evaluation of the first two terms on the right-hand side. By simple manipulations we 
deduce that 

Y{t,x) := - r a(r)/3(r,-)ip 2 (t,x,r) — 2a(r)ip(t,x,r) J /3(r,-)K(r,r)ip(t,x,r) dr^ dr 

1 / roo \ 

a{r)ip{t, x, r) I /3(r, x, r) — 2 / /3(f, -)n(r, r)ip(t, x, f) df 1 dr 


(64) 


x, r) /3(r, x,r)-2 


/3(f, -)ip(t,x,f) 


dr ) dr 


a(r)ip(t, x, r) 


ct 


P{r,-)ip(t,x,r) 


dr ) dr 


a(r)ip(t, x, r) 


d r I(t, x, r) dr. 


Next, we rewrite ip in terms of I{r) as follows: 

pt \m \ a m i,o 21{t,x,r) 

rp(r,-)ip(t,x,r) = —d r I(t,x,r) + 2r / -dr = —a r I(t, x, r) H-. 

Jr r r 

Hence, substituting this identity into (l64l) . we deduce that 


Y(t,x) = [ -^p—r/3{r,-)ip{t,x,r)d r I(t,x,r)dr 

J ro r 2 /3(r,-) 

a(r) 


r o 


(65) 


2 /3(r, •) 
a(r) 


d r I(t,x,r) ( - d r I{t,x,r) + 


2 1(t,x,r) 


dr 


OO 


r 2 /3(r, 
a(r) 


^|d r /(f,x,r)| 2 dr + f d r I 2 (t, x, r) dr 

> / */ro r rV» / 

a(r) 




r2 




/ (t,x,r 0 ). 


Thus, returning to (1631) . we substitute H(t, a;) into the first two terms on the right-hand side, and integrate by parts 
in the last term in (l()3l) recalling that r(r 0 ) = 0; this yields 

-j- [ f a(r)ip 2 (t, x, r) dr dx + 2 f f a(r)H(r)|Va;^(t,a:,r)| 2 drda; 
di Jn Jr o do Jr 0 


( 66 ) 


i(r) 


M 


/f2 Jr 0 

nOO 


r 2 P(r, 


= -2 


h' aMx ' r)] 2 irix+ 2 I,,^M 
t r »( “M 

J Q, J 


I 2 (t , x, ro) dx 


'r 0 

✓»oo 


3 /3(r, •) 


/ 2 (t, x, r) dr dx 


/ / J>(t, x)d r (r(r)a(r)) a/) 2 (t, x, r) dr dx. 

J Q J r o 


By setting a(r) := r 3 /3(r) with nonnegative /?, using Holder’s inequality for the last term, we get 


_d 

dt 


O dr 0 
+ 2 


r 3 fj(r)ip 2 (t, x, r) dr dx + 2 
r/3(r) 


r 3 /3(r)H(r)|V a ,i)’(t, x, r)| 2 dr dx 


(l Jr 0 


(67) 


/ */r 0 

/*oo 


/3(r, •) 


\d r I(t,x,r)\ 2 drdx + 2 


r2 


P{ror) 


I 2 (t, x, r 0 ) dx 


< -2 


ct 




IQ Jr 0 

+ Halloo 


/3(r, • 


/ 2 (f, x, r) dr dx 


[ r 

Q | T ( r )^( r ) , 3T ( r ) 

J £1 J ro 

* W+ «r) + r 


r 3 P(r)ip 2 (t, x, r) dr dx. 
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Consequently, if we choose p such that for all (r, it, D) £ (r-g, oo) x x M. dxd there holds 


( 68 ) 


d r 


P(r) 


> 0 , 


j3(r,u,D), 

then it follows from our assumptions on r (cf. (fT8l) l that 
d 


r ( r)d r fi(r) 


< C0), 


d t 


O J i"o 


r 6 /3(r)i/j 2 (t , x, r) dr dx + 2 


Q. J r 0 


P( r ) 

o 

r 3 p(r)A(r)\S/ x ip(t, x, r)\ 2 dr dx 


< C(K, H^olloo) / / r 3 p(r)i/j (f,x,r)drdx 

J £1 J r o 

and by Gronwall’s lemma we get 

p poo p poo 

.sup te ( 0 T ) / / r 3 p(r)ip 2 (t, x, r) dr dx + / / r 3 P(r)A(r)\S7 x i/j(t : x,r)\ 2 dr dx dt 

J Q J rn J Q J rn 


ess.: 


(69) 


<C(K, Halloo) 


O Jr 0 


r 3 f3{r)ip q (x, r) dr dx. 


In what follows, we focus on finding p that satisfies (l 68 l) . For any nondecreasing y(r) > 1 we define 

3 / \ / \ f r r]{r) dr 

p(r) := 7 (r) e- 1 ^ u ' , 

where 77 is introduced in (f20l) and we check (f68|) . For the second inequality in (j68]l . we observe that 

drP{r) dr'yir) 

-JW = — + viry ’ 

consequently, since r(r) < Cr and r] satisfies m, it suffices to choose 7 such that 

rd r y(r) 


(70) 


q(r) 


< C for all r £ [ro, 00 ] 


to ensure the validity of the second inequality in ([ 68 ) 1 . To check also the first inequality in 
aid of (l 20 l) and the fact that 7 is nondecreasing that 


d r 


P (r) 


<9r/3(r)/?(r, u , D) - /3 (r)d r /3(r, it, D) 


/3(r) 


we deduce with the 
d r (3{r) d r (3(r : u,D)\ 


/3(r,it,D) 


> 


P 2 (r, m, D) 
P(r) (d r p(r) 


P{r,u, D) \ f}(r) /3(r,it,D) 


P(r,u, D) \ ^(r) 


_ 7?(r) ) = ^(r)> 

V{ >) P(r,u, D) 7 (r) ~ 


Consequently, the first inequality in ( 1681 ) is satisfied. In addition, setting y(r) := (1 + r)“ with arbitrary a > 0, we 
see that 7 is increasing and also that m holds. Moreover, with such a choice, it is not difficult to verify using the 
definition of /? that 

7 (r) < P{r) < C* 7 (r). 

By substituting this relation into (1M1) we obtain (1521) . □ 

4. Proof of the main Theorem 

This section is devoted to the rigorous proof of our main result. To this end, we introduce several levels of 
approximation. First, since the argument is based on the energy method, we add a regularizing term to J 6 ]) as 
follows: 

d t v(t, x ) + div x (v(t, x) 0 v(t, x)) + V x q(t, x) - div x S(if>(t, x), D x v(t, x)) + e\v\ 2p '^ 2 v = f, 

di v x v(t, x) = 0 . 

Next, we also regularize the equation <[9|) for if) by adding the term —5d rr tp to its left-hand side and to improve the 
integrability we add the terms S\tp\ 2p ~ 2 tp and S\(j)\ 2p ~ 2 cj> to (J9J) and (1TTT1) respectively. Moreover, in order to prove 
the minimum principle for if) we replace i/> by its positive part 1 /’+ in several places. Also, in order to justify the 
a priori estimates for the approximating problem, we replace <f> by (</>), where the truncation function Tfc(s) is 
defined as follows: 

Tfc(s) := min(|s|, fc) sign s, s £ R. 
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Finally, we introduce > ro and consider the following equation in Q x (ro,roo), supplemented by homogeneous 
Neumann boundary conditions on (0,T) x dfl x (ro,roo) and (0,T) x f lx {ro,roo}: 


(72) 


d t if{t,x,r ) + v(t,x) ■ V x i/j{t,x,r) + T(r)T k ((/)(t,x))drip + (t,x,r) - A(r)A x ip(t,x,r) - 5d rr ip 

f'L oo 

+ 6\ili(t,x,r)\ 2p ~ 2 ip(t,x,r) = -/3(r,v,D x v)t/j{t,x,r) + 2 / /3(f, v, D x v)K(r, r)if + {t, x, f) dr. 


Similarly, we replace the integral over the semi-infinite interval (ro,oo) by one over (ro,roo) and change the corre¬ 
sponding terms in (fTT))l as follows: 

d t fi(t,x) + v(t,x) ■ V x <f>(t,x) - A 0 A x <f>(t,x ) + />|</(t,x)| 2p ~ 2 </(t,x) 

(73) f r °° f r ° f r °° 

= -T k {(j){t,x)) / d r (rT(r))ip + (t,x, r) dr + 2 / r /3(r,v,D x v)K{r,f)if + (t,x,f)drdr, 

J ro J 0 J ro 

subject to a homogeneous Neumann boundary condition on (0,T) x dfl. Having introduced these three levels of 
approximation, we then first let <5 —> 0 and r^ —> oo, and, finally, we let e —> 0. 

4.1. Existence of a solution to the e-, S-, roo-approximating problem. The existence of a solution to the 
approximating problem follows from the following result, which we state without proof, as the argument, based on 
monotone operator theory and the Aubin Lions lemma, is completely straightforward. 

Lemma 4.1. Let H C R d , d £ {2,3}, be a Lipschitz domain, r^, £ (ro,oo) and T > 0. Assume that (Al)-(A5) are 
satisfied. Moreover, let no £ L ^ div (f2), <fio £ L 2 (Ll), ipo £ L 2 (fl x (ro, Too)) and f £ L p (0, T; W}} d j p ). Then, there 
exists a triple (v,ip,(j)) such that 

(74) v £ C(0, T; L 2 (tt) d ) n L p { 0, T; W*; p iv ) n L 2p '(0, T; L 2p ' 

(75) dtv £ (i ? (0,T;^; p div )n^'(0,T;^'(O) d ))*, 

(76) </> G C(0, T; L 2 {n x (r„, rj)) n L 2 { 0, T; W 1 ’ 2 ^ x (r„, r^))) n L 2p (Q x (r 0 , rj), 

(77) d t ^£(L 2 (0,T- 1 W 1 ' 2 (nx(r 0 ,r oo )))nL 2p (Qx(r 0 ,r oo )))*, 

(78) </> £ C(0, T; L 2 (H)) n L 2 (0, T; W ll2 (fi)) n L 2p (Q), 

(79) g (L 2 (0, T; W li2 (n)) n L 2p (Q ))*, 

satisfying, for all w £ L p (0, T; W^’ div ) fl L p '(Q) d , 


(80) 


/ (d t v,w) — / n g) n • + S($(i, x), D x v(t, x)) ■ V x w + e|n| 2p 2 n-u;dxdt 

Jo Jq 

= —a* f j v-wdS+ f (f,w)dt; 

Jo Jd n Jo 


/or all w G L 2 (1F 1 ’ 2 (H x (r 0 ,roo))) n L 2p (Q x (nur*,)), 


(81) 


/ (St'i/), w) dt + / — i/jv ■ V x uj + T k ((j))Tud r ’il> + + AV x if ■ V x w + 5d r il’d r uj dr dccdi 

J 0 *^Qx(ro,roo) 

= f —5\if\ 2p ~ 2 ifu} — (3 (-,v,D x v)iJju] + 2co f fi(r,v.D x v)n{r,f)if + {f,-)dfdrdxdt-, 

jQx(r 0 ,r oo ) Jr 


for all 6 £ L 2 { 0, T; W 1>2 (fi)) n L 2p {Q), 

[ {dt<f,0)dt+ I — <fv ■ V x 6 + AoS7 x (j) ■ V x 9 + S\cl)\ 2p ^ 2 (f> 9 dx dt 
Jo Jq 

= - f T k (4>)9 f d r (rT(r))if + (r,-)drdxdt 

Jq J ro 


(82) 


and fulfilling, in addition, 


+ 2 


9r / j3(f , v, D x v)n(r, r)ip + (r, •) dr dr dec di; 


Q Jo 


v(0)=v o , cj)(0) = (j>o, 4>(0)=ip 0 
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S(0(t, x),D x v{t, x)) := I '(0(t, x), |D x v(t, a;)|) D ;E u(t, x) 


0(i,a;):= / 'y(r)i/j(t,x,r)dr. 


and 

(83) 
with 

(84) 

4.2. A priori estimates. In this subsection, we derive bounds on the solution whose existence has been stated 
in Lemma El Some of the bounds will be independent of the order of approximation; if, on the other hand, 
something depends on one of the regularization parameters, this will be clearly indicated. Moreover, since we shall 
be following, step-by-step, the formal bounds developed in Section O some of the details will be omitted for brevity. 
Also, to abbreviate the notation, we introduce f l rao := ft x (ro, Too)- We formulate the result in the following lemma. 

Lemma 4.2. Let (u,0, 0) be a solution to (1741) (|84|l constructed in Lemma \4-1\ Then, the following energy identity 
holds for all t £ (0, T): 


(85) 


^Il«( i )ll 2 + [ [ e\v\ 2p ' +S -D x vdxdr + a* [ [ |n| 2 dS'dr= f {f,v)dr+\ 
z Jo an Jo Jdn Jo 1 


Consequently, we have the following uniform a priori bound: 

pT 


(86) ess.sup tg(0jT) ||n(t)|| 2 + j ^ |M|? iP + ||S||£, + \\v\\ 2 L2{an) + e\\v\\ 2 2 p p , dt < C (J|n 0 || 2 ■ 

Moreover, if 0o and 0o are nonnegative almost everywhere, then 

(87) 0 > 0 in Q x (ro,roo), 0 > 0 in Q. 

In addition, if <fo £ then the following uniform estimate holds: 

(88) ||0||l°°(Q) < max (A' 2 , ||0 o ||oo), 
and, for all q > 3, we have the following 5-dependent bounds: 




(89) 

(90) 

(91) 


f 

Jo 


|i,2 dt <C(||0o||oo, K) 5* r- 3 + C(||0o||oo,K 


•f 

jQr 


ri/jodrdx + C(||0o||oo, K) S 2 


L 


r 3 t/jo dr dx, 


s-supt e(0 ,T) / r q ip 2 (t)drdx 

Jo. roo 

+ f f A(r)r q \V x ip\ 2 + 5r q \ip\ 2p + 5r q \d r 'il’\ 2 dr dx dt < C(q)e CST f r 9 0odrdx, 
Jo an-. 


ess.sup tg(0T) 


/ 

jQr 


r q 2 0(t)drda; < C(||0o||oo, <?, A) I S 2 r, 


(W - 6 +[ 

\ Jn, 


r q 2 0o + 5 * r 9 0 o dr da; | . 




Proof First, the energy identity (1551) directly follows by setting w := r>X[o,t] in (180[) . Then, the uniform estimate 
(1861) is a consequence of the assumption (1251) on S. To obtain the minimum principle for 0, we mimic the proof of 
Lemma IQ Thus, we set 9 := 0_ in (1821) . where 0_ := min(O,0). Hence, using the divergence-free constraint on 
the velocity v , the fact that 0+, /3 and r are nonnegative and that r is nondecreasing, we deduce that 

Therefore, as 0o is assumed to be nonnegative, we find that 0_ = 0 and 0 > 0 almost everywhere in Q. The proof 
of the inequality 0 > 0 is similar and is therefore omitted. Consequently, we can replace 0+ by 0 in (|5T1) and ([82]) . 
In order to prove ([88]) . we mimic the proof of Lemma l3~2l To this end, we begin by noting that the right-hand side 
of (1551) can be, for nonnegative 9, bounded as follows: 

<9 r (rr(r))0 + (t, x, r) dr da; dt 


[ T fc (0)d f 

J Q J r 0 

r r r o r^oo 

+ 2 9r P(r,v,D x v)n{r, r)^+(t, a;, f) dr dr dx dt 

JQ JO J rn 


<-K~ x ro / 0(T fc (0)-M) 


0(r) dr ) dx dt 
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with M := max(||^o||oo) K 2 ). Therefore, by setting 0 := (<j> — M) + in (1521) . using the fact that div^ v = 0, the 
nonnegativity of (p and the above estimate, we arrive at 

±\\^-M) + \\l<Q. 

Since cpo < M almost everywhere in the estimate (1551) directly follows. Thus, in what follows, we assume that 
k > M and therefore we can replace X). ((p) by <p in (15T1) and (15121) . 

We see that all of the above estimates are independent of the order of the approximation. In what follows, we 
shall establish several estimates that depend on some of the regularization parameters, but the estimates will become 
uniform by applying the relevant limiting procedure. Note here that while in the proofs of Lemmas 13.3113.61 we have 
formally used integration by parts with respect to r, assuming that ip vanishes sufficiently quickly at infinity, at the 
level of approximation, in a rigorous argument, such a procedure is not allowed. Fortunately, using the minimum 
principle for ip and the fact that r(r 0 ) = 0 , we can still integrate by parts in all of the desired terms at the cost of 
finally changing the equality sign to an inequality sign. 

We start with the bounds that are similar to those in Lemma 13.61 Hence, setting w := a(r)ip\[ o.sl i n (1811) with a 
suitable positive a(r), r S [r 0 , 00 ), to be selected below, we get (by recalling that H roo := ft x (r 0 , Too)) the following 
equality: 


a(r) 


(ip 2 (s) — ip^) dr dx + f ( a(r)A(r)\\7 x ip\ 2 + 8a(r)\ip\ 2p + 5a\d r ip\ 2 drdxdt 

J 0 Jn^ 


-/3(-,v,D x v)aip 2 + 2aip ( / /3(f, i>,D x v)n(r, r)ip(r) df I drdxdt 

0 J Q roo \Jr J 

— — l ( <pTad r ip 2 drdxdt — sf f a'ipd r ip dr dx dt. 

2 Jo Jci roo Jo Jci Tao 

First, we use Young’s inequality in order to absorb part of the last term into the left-hand side, to deduce the 
following inequality: 


*(r) 


(ip 2 (s) — %Pq) dr dx + j f a(r)A(r)\V x ip\ 2 + 5a(r)\ip\ 2p + ^-a\d r ip\ 2 drdxdt 

Jo Jc 2 


< 


J J —/3(-,v,D x v)aip 2 + 2aip^J /3(r,v,D x v)K(r,r)ip(r,-)drj drdxdt 


1 


4>Tad r ip 2 dr dx dt + <5 


2 a 


-ip 2 dr dx dt. 


2 Jo Jci roo Jo Jn r ^ 

Next, we integrate by parts in the third term on the right-hand side to deduce that 
/ a(r) 


Ja 

< 


(ip 2 (s) — iPq) dr dx + ( f a(r)A(r)\S7 x ip\ 2 + 5a(r)\ip\ 2p + ^-a\d r ip\ 2 dr dxdt 

Jo Jc 2 


(92) 


—y®(-, v, D x v)aip 2 + 2aip ( / /3(r,v,D x v)k(t, r)ip(f, •) df I drdxdt 

l roo r J 

- 7 ; / (p(t,x)(r(r 00 )a(r 00 )ip 2 (t,x,r 00 )-T(ro)a(ro)ip 2 (t,x,ro)) dxdt 

1 Jo J Cl 

ip 2 dr dx dt + - f f (pd r (ra)ip 2 dr dx dt 

2 Jo Jn r „ 


0 Jci r 


iq r 

2 a 


< 


J J —f3(-,v,D x v)aip 2 + 2aip /3(r,v,D x v)n(r : r)ip(r, ■) dr^j dr dx dt 
i'\ 2 1 r s r 

— ip 2 drdxdt -|— / / (pd r (ra)ip 2 dr dxdt, 

- a 2 Jo Jci T ^ 


0 JCl r 


where for the second inequality we have used the fact that r(ro) = 0 and that <p, t and a are nonnegative. We then 
follow, step-by-step, the proof of Lemma 13.61 In particular, by setting a(r) := r 3 (1 -F r) 7 e ^ r 0 ds , with 77 given by 
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4HD and nonnegative 7 , we see that the first term on the right-hand side is nonpositive and therefore 
a(r ) 


(93) 


/ 


n rc 
< S 


0 Jn 


(ip 2 (s) — ip d ) dr dx + ( f a(r)A(r)\X 7 x ip \ 2 + 5 a(r)\ip\ 2p + -a\d r ip \ 2 dr dxdt 
Jo Jn roo 

/I 2 1 r s r 

-— ip 2 dr da; dt + - / / (pd r (Ta)ip 2 dr dx dt. 

2a 2 Jo J fu 


Hence, by the assumption (TTSl) . the definition of a and Gronwall’s lemma, we see that for all q > 3 (noting that we 
set 7 := q — 3 in the definition of a) we have 


r q ip 2 (t) dr dx + 


(94) 


/0 JQ r 


r q A(r)\V x'lpl 2 + 5r q \ip\ 2p + -r q \d r ip\ 2 dr dxdt 


< C(q)e CST f r q ip d dr dx, 
J\l r 


which is nothing else than (15U1) . Next, in order to estimate the high-order moments of -0, for any q > 1 we set 


u := r q X[o,s ] hi (1511) to obtain 


(95) 


[ (r q ip(s) — r q ipo) dr dx + [ f cpTr q d r ip + Sqd r ipr q 1 dr dxdt 

Jn r „ Jo Jn r „ 


f [ m 

Jo Jn r ~. 


2p-2„/,„qr 


ipr q + /?(-, v, D x v)ipr q — 2 r q 


/3(f, v, D. r v)n(r, r)ip(f, •) df I dr dxdt. 


The first term on the right-hand side is nonpositive and can be therefore discarded. Next, we integrate by parts in 
the second term on the left-hand side and use the fact that all functions involved are nonnegative (and thus we can 
neglect their values at roo); and, finally, we use Young’s inequality in the last term on the left-hand side to deduce 
that 


(96) 


f (r q ip(s) — r q ipo) dr dx < f j C(q)S^ (6r q+2 \d r ip\ 2 + r q 4 ) + (pd r (Tr q )ip dr dx dt 

<J £2 r ^0 J Or- 


+ 


—/ 3ipr q + 2 r q 


0 Jn r 


1 0 J n ro 

P{r,v,D x v)n(r,f)ip(f,-)df\ dr dxdt. 


For the first term on the right-hand side, we use the estimate (1901) and for the remaining parts we can follow, 
verbatim, the proof of Lemma 13.51 (by using the assumption (1181) on r). Finally, we replace q by q — 2 to complete 
the proof of m- 

Finally, we focus on (l89l) . Setting 9 := (j) in (l82l) . using the bound (l 88 l) . and the assumptions (A3) and (A4), we 
deduce that 


(97) / \V X 0\ 2 + 6\0\ 2p dxdt <C(A 0 ,K, ||0 o ||oo) / / r^drdxdf 

J Q J Q Jr 0 

and (l89l) then follows by using (liiTl) with q = 3. 


□ 


4.3. The limits S —x 0+ and r^ — > 00 . This subsection is devoted to passage to the limits A —>• 0_|_ and —> 00 ; 
i.e., we eliminate the presence of the elliptic regularization of ip and we pass from the bounded interval (ro, r x ) to 
(ro,oo). We shall use the notations floo :=fix (ro,oo) and Qoo := (0, T) x Hoc. The associated existence result is 
formulated in the following lemma. 

Lemma 4.3. Let Q C R d , d £ {2,3}, be a Lipschitz domain and T > 0. Assume that (A1)-(A5) are satisfied. 
Moreover, let i>o £ L^ div , f £ L p (0,T; W~^), and let <po G L°°(H) and ipo G T^(^oo) be nonnegative; then, there 
exists a triple (v,ip,<p) such that 

(98) v £ C(0, T\ L 2 (fl) d ) n L p ( 0, T\ W^ iv ) n L 2p ' (0, T; L 2p \<A) d ), 

(99) d t v £ (L p (0,T;W^’ p d .JnL 2p \0,T;L 2p \n) d )y , 

(100) iP £ <0(0, T; L 2 (0)) n L 2 ((0, T) x (r 0 , r*); W 1 ’ 2 ^)), 

(101) dti> € (i &$1 oo) n w 1 ' 2 (n 00 ))y, 

(102) 0 £ C{ 0, T; L 2 (H)) n L 2 (0, T; W 1 ' 2 ^)) n L°°(Q), 

(103) a e 0€(L 2 (O,T;W 1 - 2 (n)))*, 
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which satisfies, for all w £ 27(0,2*; W^’j iv ) D L p ( Q) d : 

[ (d t v,w)dt— f v ® v ■ S7 x w + x), D x v(t, x)) ■ V x w + e|i;| 2p ~ 2 v ■ w cLcdt 

J o Jq 

= — a* f f vwdS+ f {f,w)dt; 

Jo JdQ. Jo 

furthermore, for all u £ L°°(0, T; fU 1,2 (f2oo) H VF 1 ’ 00 (ri 00 )) n L 2p (Q oa ) fulfilling w(oo, •) = 0 one has: 


(104) 


(105) 


/ (dtf-’, co) dt + / — ipv ■ \7 x uj — 09 r (ru;)'0 + A'Vxif • V^wdrdxdt 

Jo J Qoo 

= J —/3(-,v,D x v)fi>oj + 2uj J /3(r,v,D x v)K,(r,f)'ip(r,-)drdrdxdt; 


and, for all 9 £ L 2 {0,T;W 1 ’ 2 (fl)) one has: 


(106) 


f (dt<j>, 9) dt + I —fiv-V x 9 + AqS 7 x (f-V x 9 dxdt 
Jo Jq 

= — f 4>6 f d r (rT(r))il’(r, •) dr dxdt 

Jq Jr 0 

r r^o roo 

+ 2 9r /3(r,v,D x v)n(r, r)if(r, •) dr dr dx dt; 

Jq Jo Jrn 


’ Q JO J r 0 

with v, cj) and if attaining the initial data v 0 , cj) 0 and ipo, respectively, in the sense that 
(107) limsup || v(t) - uolU + || (j>{t) - fio\\l + || ip(t) - ifoWl = 0, 

t —^0+ 


and fulfilling (1831) . Moreover, the solution satisfies ©, ®. (EZD and (1881) . In addition, for all q> 3 we have the 
following uniform estimates: 

(108) f ||0||i,2di < C'dl^olloo.-K’) [ ripodrdxdt, 

JO J Qoo 

(109) ess.sup tg ( 0T ) f r q ip 2 (t) dr dx + ( [ A(r)r q \V x ip\ 2 dr dxdt < C{q) f r^odrdx, 

' Jn oo Jo JJ sioo 

(110) ess.sup te i 0 , T ) f r q ~ 2 if{t)drdx<C(\\(j) 0 \\ oo ,q,K)f r q ~ 2 if 0 dr dx. 

J^oo JQoo 

Proof. To prove Lemma 14.31 we use the existence result obtained in Lemma 14.11 in conjunction with the uniform 
estimates derived in Lemma 14.21 Thus we set <5 := n^ 1 and Too := Inn in Lemma 14.II and denote the corresponding 
solution by (v n , , tp n ). Using the estimates (1511 (15T1) . the fact that ipo is a smooth compactly supported function, 
and defining all functions involved to be identically zero outside (ro,^), we deduce the existence of a subsequence 
that we do not relabel and the existence of (v, <f>, if, S), where (f> and if) are nonnegative functions, such that 


(111) 

v n V 

weakly in 27(0, T; W^ p div ) n L 2p {Q ) d , 

(112) 

S n ->• S 

weakly in L p ( Q ) dxd , 

(113) 

(j) n -»* 0 

weakly* in L°°(Q), 

(114) 


weakly in 27(0, T; kL 1,2 (U)), 

(115) 

l/) n 

weakly in L 2 (<5oo) nl 2 (0,T; L 2 oc (r 0 , oo; fU 1,2 (U))) 
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Moreover, using the weak lower semicontinuity of the norm function, and the special choices of S and made 
above, we deduce that the following uniform estimates hold for q > 3: 

(116) Ml~(Q) <C(K, \\Moo), 

(117) ess.sup tg(0jT) ||u(f)|| 2 + f |Mli iP + ||S||p/ + |Mli, 2 (afl) + e||u||pf dt < C(vq, /), 

Jo 

(118) [ ||0|| 1,2 dt < (7(11001|oo, it") f rt/jQ dr da: dt, 

Jo Jn a o 

(119) ess.sup t g( 0 ,t) [ r 9 0 2 (t) dr da; + f f A(r)r q \V x ijj\ 2 dr dx dt < C (q) f r 9 0gdrda;, 

Jo ,oo Jo Jq , 0 o ^oo 

(120) ess.sup tg ( 0 ,T) [ r q - 2 i/j(t) dr dx < Cdl^oWoo, q, K) f r q ~ 2 i/j 0 dr dx. 

Our objective is to let n —> oo (and consequently <5 —>■ 0+ and — > oo)) in (|80l) (f82]> . To do so, we first observe 

that (l80l) (l82)l imply the following weak convergence results: 


(121) 

d t v n d t v 

weakly in 

(LP(0,T-,W^ div )nL 2 v(Q) d y, 

(122) 

5 t 0" - d t 0 

weakly in 

(i 2 (0,T; H rl ’ 2 (fl)))*, 

(123) 

9 t 0” ->■ dti> 

weakly in 

(L 9 (0, T; W 1,2 (K) D W 1,q (K)))* 


for sufficiently large q and any K := x (ro, k) with arbitrary fixed fcgR. In addition, red by noting the uniform 
(with respect to r^) bounds (lll(il) (11201) . we see that the weak limit dtip satisfies (11011) . Consequently, we can use 
the Aubin-Lions lemma to deduce strong convergence, and then extract subsequences still labelled by the index n 
(i.e., without indicating the subsequences in our notation), such that 

(124) v n —t v a.e. in Q, 

(125) 0 71 —y 0 a.e. in Q. 


z(r)%l>(t , x, r) dr weakly in L 2 (Q). 


However, we cannot claim the same convergence result for 0" because of the lack of the compactness with respect 
to r. Nevertheless, we will show that, for arbitrary z G 2?(ro, oo), we have 
/■rS, r<x> 

(126) / z{r)ip n {t, x, r) dr —> / z(r)ijj(t,x,r) dr strongly in L 2 (Q). 

Jr 0 Jr 0 

First of all, it follows from (11151) that 

pr n 

(127) / z(r)i/j n (t, x, r) dr - 
J r o 

Hence the limit is defined uniquely. Next, denoting 

(128) <ip*(t,x):=[ z(r)ip n (t, x,r) dr, 

J r 0 

we can set in m ui := z(r)ip(t,x), where ip G T>(Cl) is arbitrary, and for sufficiently large n such that In n is not in 
the support of z we obtain the identity 

J (dti^g, <P^ dt + J • S7xW + V x 0^ z ■ X7 x tp da; dt 

(129) = J <P—S\ip n \ 2p ~ 2 ip n z — /3(-, v n , D x v n )ip n + 2z /3(f, v n , D x v n ) K (r, r)0 n (r, •) df) drj dxdt 
— J (f)ip —d r (Tz)i/) n dr'j+Sip^J dr^drzdr^j da;dt. 

Consequently, using the a priori estimates (1M1) (1M|) . we see that, for sufficiently large q, 

(130) —*■ dt'ipz weakly in T 9 (0, T; W 0 1,g (f])). 

On the other hand, using (li!Ul) we also have that 


||C||? 2 d t<C(z) [ |0"| 2 + A(r)r 3 |V x 0"| 2 drdxdt < C, 

Jq 
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where the first inequality follows from the fact that 2 has a compact support. Hence the Aubin-Lions lemma 
completes the proof of (11261) . 

Next, we shall apply a similar convergence argument to tl> n 1 which is defined in (| 8 |). Since 7 is a continuous 
function, we can find a sequence of G 2d(ro, 00 ) such that 7 E /* 7 almost everywhere and also locally in C(r 0 , 00 ). 
For such an approximation, we can however use the convergence result (11261) and obtain 

f r °° f°° 

(131) / 'y s (r)ip n (t,x, r) dr —>• / 7 E (r)^>(t, x, r) dr strongly in L 2 (Q). 

J r 0 J r 0 

Hence, for the original function J> n we have that 

f dxdt := f f 'y(r)(tp n (r, t, x) — ip(r, t, x)) dr 

JQ Jq Jr n 


dx dt 


(132) 


< 2 


( 7 ( r ) _ 7 e (r))(V>”(r, t, x) - t, x)) dr 


dx dt 


+ 2 [ [ 7 e (r)(V’"'(r,t,x) --0(r,t,x))dr 
Jq Jr 0 


dx dt. 


Consequently, using (11311) . we deduce that 

lim / \J) n — xj)\ 2 dxdt < 2 lim [ [ (y(r) — 7 E (r))(V>"(r, t, x) — ijj(r, t, x)) dr 

' n^-oo Jq J 


(133) 


< 4 lim 

n—>■ 00 


( 7 ( r ) - 7 e (r)) (J) n (r,t,x) -• 0 (r,t,x))dr 


dx dt 


dx dt 


+ 4 lim / [ | 7 (r)||^ n (r, t, x) — ^(r, t, x)| dr 

n ~^°° JQ J r . 


dx dt. 


Thanks to the uniform convergence of y E on compact sets, we can also easily let e -A 0+ (with the aid of the uniform 
bound (l9Tl) l to get 

p 00 2 


n—too 


n—t 00 


(134) 


lim / \ip n — J)\ 2 dxdt < 4 lim 
'Q 

poo 

< 


Q 


| 7 (r)||V>"(r,t,x) -ip(r,t,x)\dr 

dr r r°° 

4 / —— dr lim / / | 7 (r)| 2 r 2 |i/’”(r, t,x) — i/j(r,t,x)\ 2 drdxdt 

Jr * r n ^°° Jo Jr * 


dx dt 


< K(r*) 1 lim [ [ r 2e+2 (\ip n (r,t,x)\ 2 + \4>{r,t,x)\ 2 ) dr dxdt < C(K,i/jo){r*) 1 , 

n - > °° do dr* 


1 Q Jr 

where we have used (I231l . the fact that if>o is compactly supported and the estimate (19U1) . Consequently, letting 
r* -> 00 we deduce that 

(135) J) strongly in L 2 (Q). 

We now have all ingredients in place to complete the proof. First, it is standard to let n —> 00 in (l80l) to deduce, 
for all w G U>( 0, T; W^ iv ) n L 2p (Q) d , that 

/ (d t v,w) — / v ® v ■ V x w + S • X7 x w + e\v\ 2p ~ 2 v ■ w dxdt 
do Jq 

= — a* f f v-wdS+ j (f,w)dt, 
do Jon do 

as well as the first limit in (I107D . Consequently, setting w := v, we obtain the following identity: 

(137) \b{T)\\l+ j f e\v\ 2p '+S-D x vdxdr + a* f [ \v\ 2 dSdr = f (/, v) dr + ^||v 0 |||. 

z do dn do dan do z 

Then, using the weak lower semicontinuity of the norm function, the energy identity (1851) with t = T and the weak 
convergence result (11111) . we deduce that 


(136) 


(138) 


limsup / S(Jj n ,D x v n ) ■ D x v n dxdt < / S D^vdxdt. 
n-y 00 Jq Jq 
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Moreover, using (11351) . Lebesgue’s dominated convergence theorem and the assumption (f2o|) . we also see that 
(139) S{Jp n ,D x v)^S{iP,D x v) strongly in L p '(Q) dxd , 

which, when combined with (11381) and (11111) , (11121) , leads to 


(140) 


limsup f (S^ 71 ^xV 71 ) — S(ijj n ,D x v)) ■ (D x v n — D x v) dxdt = 0. 

n —loo J Q 


The strict monotonicity assumption (1251) then implies that there is a subsequence such that 


(141) D x v n —>■ D x v a.e. in Q, 

and consequently we have that S = S(jp,D x v). Finally, in view of the convergence results and a priori estimates 
obtained, one can let n —> oo in m, m to deduce (11051) . (11061) . We note in this respect that, thanks to the 
nonnegativity of both %p (and cp) the + symbol, indicating the nonnegative part of a function, can be omitted from 
ED- Furthermore (1811) is strongly nonlinear because of the presence of the term 5\ip\ 2p 2 tp, but we are considering 
the limit S —> 0 + , so this term vanishes in the limit of <5 —> 0 + thanks to the a priori estimate (IDT71) . 

Finally, by using a standard parabolic regularity result, one can show the attainment of the initial datum (11071) 
for cp, and also that ip G C w (0, T; L 2 (f2 oc )) fulfils, for t —> 0+, 

(142) ip(t) —^ ip o weakly in L 2 (f2oo). 

To strengthen this convergence result, we recall (ED, which, thanks to Gronwall’s lemma, leads to 

( 143 ) f a(r)(ip n (s)) 2 dr dx < e c( ~i^o, K )t f a{r)ipQ dr dx. 

J n rSb 

Hence, letting n —> oo and using the weak lower semicontinuity of the norm function, we get 

f a(r)(ip(s)) 2 dr dx < f a(r)ipQ dr dx. 

Jflaa Jfloa 


This directly leads to 


limsup / a(r)(ip(s)) 2 dr dx < 

>0+ JO. oo 


/ a(r)ipQ dr dx, 

J Qqq 


which, when combined with (11421) . yields (11071) for ip. 


□ 


4.4. The limit e —> 0+. In this final subsection, we complete the proof of the main theorem in the paper. For this 
purpose, we use the existence result from Lemma l4~3l Hence, for tpo g T 1 (H; L\» (r 0 , oo)) D T 2 (H; Lg. (r 0 , oo)) with 
0* > 9 > 1 and 0\ > 3 we find a sequence tpg G X^floo) that converges strongly to ipo in the corresponding spaces. 
We then denote by (v e ,(p s ,ip e ) the solution constructed in Lemma l4~3l with the initial data (po and ip^. Our goal 
is now to let e —> 0+ and obtain a solution whose existence is claimed in Theorem 11.11 Henceforth, we denote by 


C a generic constant that may depend only on the data 
previous section, we have that (1571) is valid and 


(144) 

ess.sup te(0T) 

\\v E (t)\\ 2 + [ ||n e ||? !p + 
Jo 

(145) 

\W\\l°°(Q) < 

c, 

(146) 

ess.sup t6(0T) 

[ r e1 > (ip e ) 2 (t) dr dx + 
Jo^ 

(147) 

ess.sup te(0T) 

j r Bl ip £ (t) dr dec < C. 

J ^oo 


but not on e. Recalling the estimates established in the 

\\S E \\i + \\v s \\l 2 {dn) +e\\v^dt<C, 

[ [A(r)r e ^\V x tp s \ 2 dr dx dt < C, 

I o Jo^ 
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Therefore, using the same arguments as before, one can deduce that there exists a quadruple (v, S, i/>, <j>) such that, 
for sufficiently large q > 1 , one has 


(148) 

v E 

—^* V 

weakly* in L p ( 0, T; W^ iv ) n T°°(0, T; L 2 (n) d ), 

(149) 

S E 

-S 

weakly in L p '(Q) dxd , 

(150) 


-*<t> 

weakly* in L°°{Q), 

(151) 



weakly in L 2 (0, T; W l,2 (Vi)), 

(152) 

r 


weakly in L 2 (Q 00 ) DL 2 ( 0 ,T; if oc (r 0 , 00 ; lT 1 ’ 2 (fl))), 

(153) 

r 


weakly* in L°°(0,T; L\.(r Q , 00 ))) D i°°(0, T; L 2 (fl; Lg. (r 0 , 00 ))), 

(154) 

v E 

—*■ V 

weakly in L 2 (0, T; L 2 (dfl) d ), 

(155) 

d t v E 

d t v 

weakly in L 9 (0, T; W“^), 

(156) 

d t f 

dt<\> 

weakly in (L 2 (0, T; W 1 ’ 2 ^)))*, 

(157) 

dt*l) E 

dti> 

weakly in (L ? (0, T; W 1 , 2 (r2<x>) D kF 1 , 5 f (fl 00 )))*, 

(158) 

v n 

-A V 

a.e. in Q , 

(159) 

r 


a.e. in Q. 


Based on these convergence results, we can deduce similarly as before that, for a sequence of jg £ V(ro,oo) such 
that 7 g 7 , almost everywhere and also locally in C(ro,oo), we have 


(160) 


poo POO 

/ 'ys{r)t/) E (t, x, r) dr —> / x, r) dr strongly in Lr(Q). 

J rn J rn 


Next, we slightly change the convergence result for the original sequence since the initial data are assumed to be 
only in Lj.. Thus, we focus only on L 1 convergence; in particular, we note that 


Jq 


— ibldxdt := 


l{r){ip £ (r, t, x) — t, x)) dr 


(161) 


< 


' r 0 
POO 


dx d t 


( 7 ( r ) - 75 (r))('i/’ £ (r, t, x) - i/){r , t, x)) dr 


'r 0 


dx d£ 


7 g(r)(ip e (r, t, x) - t, x)) dr 


dx df. 


Consequently, using (11601) . we see that 


lim / \i/j e — i/j\ dx dt < lim 


(162) 


£—»■() + 


’Q 


£—^0-|- ^ 


+ lim 

e ^o+. 


( 7 ( r ) _ 7 s (r))(ip E (r, t, x) - ip(r, t, x)) dr 


Q 


p OO 

/ \l{r)\\ij} £ {r,t,x) - ip(r,t,x)\dr 

J r* 


dx df 


dx d t, 


and thanks to the uniform convergence of 7 g on compact sets, we can also easily let <5 —> 0+ (with the help of the 
uniform bound (11471) 1 to deduce using (HM1) that 


(163) 


lim / \ij) E — ij)\dxdt < lim / / | 7 (r)||' 0 £ ( r , t, x) — ip(r, t, x)| dr 

+Jq e ^°+J Q Jr * 

p OO 

/ r 9 \ip E (r, t, x) — t, x)| dr 
J r* 


dx d t 


< C lim 
£ —^ 0 -|- 


dx d t 


< C(r*) e - e H\m\Ll, + \m\Ll.) < C(r *) e ~ e ' r *l°° 0, 


and consequently 
(164) 

Based on these convergence results, it is now an easy task to let e —> 0+ in (11041) (11061) to get 
that we can show that 


ij) e -A if) strongly in ^(Q). 


provided 


(165) 


D x .n e -A D x v strongly in L 1 (Q) dxd . 
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In addition, the proof of the attainment of the initial data for v and (f> is rather standard, and for i/j we can use the 
same scheme as in the previous section. Also to prove the conservation of mass identity for the polymer chains, we 
follow the proof of Lemma 13.31 but with a proper cut-off function. Specifically, we can set 

z ■= X[o,s] in m and P '■= X[o ,s]fr]k(r) in ([331), 

with arbitrary rj k £ 2?(K), and after summing the resulting identities we deduce that, for almost all t £ (0,f), we 
have 

/ J>(s,x) - </> 0 (x) + I / rrj k (ip(s,x,r) - r))dr ) dx 
J Q \J vq / 

= — / / r77/c^(r, x, r) dxdr dt 

Jo Jq^ 

n /*oo 

rrj k / n(r,r)f3(f,-)il)(t,x,r)drdxdrdt 
loo Jr 

/ / (d r (rTTih) — d r (rT))<j>(t, x)ip(t, x, r) dxdr dt 

Jo J 

/»00 

r / /3(r, -)n{r, x, f) df drdxdf. 


(166) 


o Jn Jo 


Next, we assume that ^ e 1 on (0,ro) and we follow the proof of Lemma [3.31 to get, with the help of (l22l) . the 
following chain of equalities: 

n f> OO 

rrjk / n(r,r)l3(r,’)ip(t,x,r)drdxdrdt 

} -oo Jr 

m r 0 poo 

r / j3(f,-)K(r,f)i(>(t,x,f)drdrdxdt 
Jr 0 

n nOO 

rr] k / n(r,r)(3(r,-)i/j(t,x,f)dfdxdrdt 

loo J r 

m VQ poo 

rr lk / /3(r, r)ip(t, x, f) dr dr dx df 

J r 

m OO n OO 

rrj k (r) / K(r,r)/3(r,-)i/}(t,x,r)drdxdrdt 

Jr 

/ rrjk(r)n(r,f)dr \ p(f,-)Jj(t,x,f)dfdxdt 


= 2 


= 2 


= 2 


o Jn Jo 


= 2 


OO / poo 


0 J Q J ro \J 0 


rrik{r)n{r, f) drj /3(f, -)^( t , x, f) df dx d< 

m oo 

fi?fc(f)/3(f, ')i/)(t,x,r) df dx dt 

0 

J r 2 rj k (r)drj /3(r,-)ijj(t,x,f)dfdxdt. 


0 J f2 J ro 


By substituting this identity into (11661) . we obtain 

J <j>(s,x) - <f>o(x) + {^J rr) k {i(>{s,x,r) - ip 0 {x,r)) dr^ dx 
/ / (rrr] k (r) + (r] k (r) — l)d r (rT))(j)(t,x)'il>(t,x,r)dxdr dt 

Jo Jn a o 

(J r 2 rj k (r)dr\ ^(f,-)ip{t,x,f)drdxdt. 


(167) 


0 J 17 J ro 


Next, we let rj k /■ 1. Thus, we set t] k {r) = 1 for r < k and rj k (r) = 0 for r > 2k such that \r/' k \ < Ck 1 . Using 
this definition it is not difficult to let k —>• oo in the terms on the left-hand side of (11671) by using the monotone 
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convergence theorem. For the term on the right-hand side, we can use (TTH1) to deduce that 

n (rrri' k (r ) + ( r] k (r ) - l)d r {rT))4>(t,x)tp(t,x,r) dxdr dt 

!„ 

f~ l IJ r 2 r]' k (r)dr\ (3(f,-) r ip(t,x,f)drdxdt 


0 J Cl J r*o 


< 


m oo r s p poo 

</>(f, x)ip(t, x, r) dxdr dt + C / / / fip(t,x,r)dfdxdt 

Jo J Cl J k 


< C(IIV’lloo) / / / ftp(t,x,f)dfdxdt k ^°0, 

Jo J Cl J k 

where the last convergence follows from the fact that ip £ T)(<3oo)- 

Hence, it remains to show (116511 . First, following (3.56), (3.57) and (3.60) in |8], using the fact that H £ C 1 ’ 1 we 
can find q\ : gf and some q* > 1 such that 

(168) q\ —qi weakly in L p (Q), 

(169) 92^92 weakly in L q *{Q), 

(170) gf — > 52 strongly in L h (Q) for all h £ [l,g*), 

fulfilling, for all w £ L°°(0, T; H /1 ’ 00 (f2) d fl W^’ 1 ), 

[ ( d t v e ,w)— I v E ® v E ■ S7 x w + S('0 £ , D x v e ) ■ S7 x w + e\v e \ 2p ~ 2 v e ■ w dxdt 
Jo Jo 


(171) 


[ f v e ■ w dS + f {f,w)dt+ f (gf + gf) div x wdxdt 
Jo J dCl Jo J Q 


and using the convergence results (114811 (115911 , we also get 


(172) 


(d t v,w) — / v ® v ■ y x w + S • S7 x w dx dt 

Jq 


I f v-wdS+ f (f,w)dt+ f (gi + g 2 ) div^, w dxdt 
Jo Jdn Jo Jq 


which is nothing else than (1521) with q given as q := q\ + g 2 provided we show (116511 to identify S. We now set 
n := [e -1 ], reinstate the index n for all functions concerned, and define 


f n :=- \v n \ 2p -V, 


G” 

H" 

H 


= v n <8> v n - v ® v + (gf — g 2 )l, 
= -S(^ n ,D I «") + gl*l, 

= -S + gil. 


Hence, we may apply Lemma l2Jl Thus, we set A* := k and we see that we can find a sequence X k £ [k,Ck pk ] 
and the corresponding sequence u n,k fulfilling (15U1) . (TUI) . Next, for any nonnegative g £ T>(Q) we set w := u n,k in 
(I171D . (1 172 [1 : by subtracting the resulting equations, using (151)1) . (1571) and (|55|) . we obtain 

(173) limsup f (S(t/>", Djw") — S) • D x (u nk )g — g(q™ — gi) div^, u n ' k dxdt < C(g)(k 1 ~ p + fc -1 )' 3 . 

n->oo J Q 

Moreover, by using (SQl), (SB), and also the fact that 

S$ n ,D x v) ->S$,D x v) strongly in L p (Q ) dxd , 

which is the consequence of (11641) . Lebesgue’s dominated convergence theorem and the assumption (1551) . we see 
that (11731) reduces to 

(174) limsup / (S(^> n , b x v n ) — S(V>", D x n)) • D x (u n,k )g — gg" div^, u n,k dxdt < C{g)(k 1 ~ p + fc -1 )' 3 , 

n—>• oo J Q 
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and consequently, after denoting by Q g the support of g and using the definition of E%, we have that 

limsup f (S(ijj n ,D x v n )—S('ip n ,D x v)) ■ D x (v n — v)gdxdt 
n—>oo J Q g \E£ 

O 75 ) <C(g) limsup/ (|H"| + |H|)|D a; (tt ri '’ fe )| dzdf + C(g)(k 1 ~ p + k~ l Y 

n—too JQ g CE™ 

< C(g)(k 1 ~ p + k~P) + C{g)(k 1 ~ p + k~ x f. 

Finally, using the monotonicity of S, Holder’s inequality and (l4dl) . we have that 

limsup f \J (S 0 x v n ) — SD a ,t>)) • 0 x (y n — v)gdxdt 

n—too J Q 

< lim sup I J (S(^ n , D x v n ) — S('0", D^u)) ■ — v)g da: dt 

n->oo J Q g \E£ 

(176) + lim sup I \J (S^jDjU™) - S(ip n , D. r i>)) • D x (v n - v)gdxdt 

n->oc jQ g nE£ 

< C lim sup ( / (S(i/j n ,D x v n )-S(ip n ,D x v)) ■ D x (v n - u^dxdf) +C\Q g ^E%\^ 

n ^°° \ J Q a \Ek J 

< C(g)(k 1 ~ p + k~ 0 )? + C(g)(k 1 ~ p + k~^ + Ck~% k ^° 0. 

Thus, using the strict monotonicity of S, see ([25]) . we deduce (for a subsequence) (|165|h which completes the proof. 

□ 


Remark. In order to keep the length of the paper within reason, we focused here on the mathematical analysis of 
a model that only admits polymerization between a polymer and a monomer. The mathematical analysis of a more 
complicated model, which also includes polymerization between different polymer chains, is of interest. It would 
involve coupling the generalized Navier-Stokes system of a form considered here, with a continuous coagulation- 
fragmentation equation, where the viscosity coefficient in the Navier-Stokes momentum equation depends on the 
size distribution function satisfying a coagulation-fragmentation equation whose nonlocal terms account for the 
formation of polymer chains by coalescence of smaller chains, the breakage of polymer chains into two smaller 
pieces, the depletion of polymer chains by coagulation with other polymer chains, and the gain of polymer chains as 
a result of the fragmentation of larger chains. We expect though that the analysis of such a more complicated model 
would proceed along similar lines. 
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